Hydrodynamics of the double-wave structure of insect spermatozoa flagella 
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In addition to conventional planar and helical flagellar waves, insect sperm flagella have also been 
observed to display a double-wave structure characterized by the presence of two superimposed 
helical waves. In this paper, we present a hydrodynamic investigation of the locomotion of insect 
spermatozoa exhibiting the double-wave structure, idealized here as superhelical waves. Resolving 
the hydrodynamic interactions with a non-local slender body theory, we predict the swimming kine- 
matics of these superhelical swimmers based on experimentally collected geometric and kinematic 
data. Our consideration provides insight into the relative contributions of the major and minor 
helical waves to swimming; namely, propulsion is due primarily to the minor wave, with negligible 
contribution from the major wave. We also explore the dependence of the propulsion speed on 
geometric and kinematic parameters, revealing counter-intuitive results, particularly for the case 
when the minor and major helical structures are of opposite chirality. 



I. INTRODUCTION 

Locomotion in fluids is ubiquitous in nature, with ex- 
amples spanning a wide range in size from bacterial 
motility to the swimming of whales. It plays fundamental 
roles throughout the lives of animals in such endeavors as 
prcdation and finding a mate for reproduction [1]. Lying 
along the interface between biology and fluid dynamics, 
biological locomotion at small scales has received sub- 
stantial attention from biologists and engineers in recent 
years [2, 3]. 

The physics governing locomotion in fluids is very dif- 
ferent for microscopic organisms (e.g. bacteria, sperma- 
tozoa) and macroscopic organisms (e.g. fish, humans). 
The dramatic difference is due to the competition be- 
tween inertial and viscous effects in the fluid medium. 
The Reynolds number, Re = U L/v, (with U and L char- 
acteristic velocity and length scales, and v the kinematic 
viscosity) is a dimensionlcss parameter which measures 
the relative importance of the inertial forces to viscous 
forces in a fluid. Locomotion of larger animals in flu- 
ids takes place at moderate to large Reynolds numbers, 
where inertial forces dominate. At this scale, swimming 
and flying are generally accomplished by imparting mo- 
mentum into the fluid opposite the direction of locomo- 
tion. Microorganisms meanwhile inhabit in a world of 
low Reynolds numbers, where inertia plays a negligible 
role and viscous damping is paramount. The Reynolds 



number ranges from 10 -6 for bacteria to 10 -2 for sper- 
matozoa [4]. The absence of inertia imposes stringent 
constraints on a microorganism's locomotive capabilities. 

Many microorganisms propel themselves by propagat- 
ing travelling waves along one or many slender flagella 
[4] . The motility features of these flagella depend on the 
cell type, either prokaryotic (cells without a nucleus) or 
cukaryotic (cells with nuclei) . The flagella of prokaryotic 
bacteria, such as those utilized by Esherichia coli, are he- 
lical in shape and are passively rotated at their base by 
a motor embedded in the cell wall. The rotation propa- 
gates an apparent helical wave from the sperm head to 
the distal end of the flagellum, propelling the cell in the 
opposite direction. Eukaryotic flagella exhibit a differ- 
ent internal structure, called an axoneme, which is com- 
posed of microtubules, proteins and protein complexes 
such as dynein molecular motors. The dynein arms con- 
vert chemical energy contained in ATP into mechanical 
energy, inducing active relative sliding between the mi- 
crotubules, which in turn leads to bending deformations 
which propagate along the flagellum. A common struc- 
ture of the axoneme has a ring of nine microtubule dou- 
blets spaced around the circumference and two additional 
central microtubules (the so-called 9+2 axoneme). Other 
variations of the axonemal structure have also been ob- 
served [5]. 

Generally, three levels of complexity in undulatory 
beat patterns are observed in eukaryotic flagella [5, 7], 
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FIG. 1: A hierarchy of the complexity of flagellar beating pat- 
tern observed in eukaryotic cells, (a) Planar-wave pattern in sea- 
urchin spermatozoa flagella [6]; (b) helical- wave pattern in Gryl- 
lotalpa gryllotalpa [7]; (c) double- wave pattern in Haematopinus 
suis [7]. All images reproduced with permission; (a) from R. Rik- 
menspocl and C. A. Isles, Biophys. J., 47, 395-410, 1985, copyright 
1985 Elsevier; (b) & (c) from B. Baccctti, Adv. Insect Physiol., 9, 
315-397, 1972, copyright 1972 Elsevier. 



following a hierarchy in the structure of the axoneme: 

(1) the lowest in the hierarchy is a simple planar beating 
pattern, as in human and sea-urchin spermatozoa flag- 
ella (Fig. la), with the common 9+2 axoneme structure; 

(2) a more complicated three-dimensional helical beating 
pattern is exhibited by some insect spermatozoa with a 
9+9+2 axoneme, as in Gryllotalpa gryllotalpa (Fig. lb); 

(3) the highest level of complexity is a double-wave pat- 
tern observed in some insect spermatozoa with a 9+9+2 
axoneme and accessory bodies endowed with ATPase ac- 
tivity, as in Haematopinus suis (Fig. lc). A vast diver- 
sity in sperm structure is found in insects [8], and the 



hierarchy described is also observed even just within the 
realm of insect spermatozoa flagella [5, 7]. Fig. 2 shows 
a schematic diagram of an pterygote insect flagellospcrm 
and its ultrastructure (reproduced with permission from 
Ref. [5]). Similar to the 9+2 axoneme observed in cilia 
and flagella of many plant and animal cells, the central 
core of the insect sperm axoneme is composed of two 
central microtubules surrounded by a ring of nine micro- 
tubule doublets. However, the ring of nine microtubule 
doublets is surrounded by another nine accessory tubules, 
forming the characteristic 9+9+2 arrangement of the in- 
sect sperm axoneme. In addition to the more complicated 
microtubule arrangement, two prominent features of in- 
set spermatozoa flagella are the mitochondrial derivatives 
and accessory bodies running along the axoneme (see 
Ref. [5] for a thorough review of insect sperm structure) . 

Although the structure of many different spermatozoa 
has been examined, the rapid and divergent evolution in 
sperm morphology is not well understood [5, 9]. Hydro- 
dynamic considerations of the relationship between flag- 
ellar morphology and functional parameters such as the 
swimming speed may provide useful information for ex- 
plaining the evolutionary divergence. Due to its intricate 
nature, the double-wave structure is less well explored 
than wave types sitting lower in the hierarchy. Rele- 
vant studies on the planar and helical wave structures 
are abundant and well developed (see the classical and 
recent reviews [2-4] ) , but we are not aware of any hydro- 
dynamic studies dedicated to the double- wave structure. 
Here we present a hydrodynamic study on the motility 
of insect spermatozoa exhibiting a double- wave beat pat- 
tern. 

The double-wave pattern is characterized by the si- 
multaneous presence of two kinds of waves, a minor wave 
with small amplitude and high frequency superimposed 
on a major helical wave of large amplitude and low fre- 
quency. The minor wave has also been observed to be 
approximately helical [10], and the combined activity 
of the two is described as a double-helical beating pat- 
tern [5]. The double- wave structure was first observed 
in Tenebrio molitor and Bacillus rossius by Baccetti et 
al. [11, 12], and was also later found in Lygaeus [13], 
Culicoides melleus [14], Aedes notoscriptus [15], Cer- 
atitis capitata [16], Drosophila obscura [17], Megaselia 
scalaris [18], and more recently in Aleochara curtula [10] 
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FIG. 2: Schematic representation of the ground plan of an pterygote insect flagellosperm and its ultrastructurc. (A) A typical filiform 
insect spermatozoon. Although not easily visible from the outside it can be divided into five distinct parts: acrosome, nucleus, centriole 
adjunct, flagellum, and tail end. (B) Cross section of the acrosome showing its trilayered arrangement of an inner acrosomal rod (ar), 
an acrosomal vesicle (av) and an outer extra acrosomal layer (el). (C) Cross section of the nucleus (n) showing condensed chromatin. 

(D) Cross section of the posterior centriole adjunct region. This part of the spermatozoon is characterized by the electron dense centriole 
adjunct material (ca), often surrounding the anterior part of the axoneme (ax) and the tip of one of the mitochondrial derivatives (md). 

(E) Cross section through a representative segment of the flagellum. In addition to the axoneme, two accessory bodies (ab) and the 
two mitochondrial derivatives of often different size (mdl, md2) can be seen. The mitochondrial derivatives typically bear paracrytallinc 
inclusions (pc). (F) Cross section through the tail end showing dissociated axonemal tubules. (G) Cross sectional representation of a 
typical 9] 9] 2 insect sperm axoneme. nine microtubule doublets (d) with associated dynein arms (da) and radial spokes (rs) are connected 
to two central microtubules (ct) via the central sheath (cs). The doublets are in turn surrounded by nine accessory tubules (at). Accessory 
tubules and doublets are linked together by intertubular material (itm). (H) Schematic cross sectional drawing of an axonemal doublet 
showing the protofilament arrangement of the A and B subtubules. The radial spoke (rs), the inner dynein arm (ia) and the outer dynein 
arm (oa) are attached to the A subtubule. This figure and caption are reproduced with permission from M. Werner and L. W. Simmons, 
Biol. Rev., 83, 191-208, 2008, copyright 2008 John Wiley and Sons. 



and Drusilla canaliculata [19]. Figure 3 compiles a col- 
lection of images of spermatozoa exhibiting the double- 
wave structure. Werner & Simmons [5] have presented a 
thorough review of this complex structure in insect sper- 
matozoa. 

Most studies on insect spermatozoa focus on the sperm 
ultrastructure, and there are very few studies on insect 
sperm motility [5]. Many important geometric and kine- 
matic data required for hydrodynamic modeling of the 
double- wave structure are unavailable. In particular, we 
are not aware of any information about the chirality of 
the minor helical structure relative to the major helical 
structure. The generation and propagation mechanism of 



the double- wave is also not yet fully understood. Baccetti 
et al. [11, 12] have suggested that the accessory bodies 
and the axoneme are responsible for the major and minor 
waves, respectively, whereas Swan [15] has stated that the 
major wave may be due to the sliding of the accessory 
tubule against the axonemal doublets. More recently, 
Werner et al. [10] have proposed a completely different 
line of thought, suggesting that the major wave is not 
in fact a real wave but a static helical structure formed 
due to the coupling of static forces of the axoneme, mi- 
tochondrial derivatives, and plasma membrane. The ap- 
parent propagation of the major wave could be due to 
the passive rolling of the entire cell and might in fact be 




FIG. 3: The double- wave beating pattern observed in insect spermatozoa fiagella of different species: (a) Megaselia scalaris [18]; (b) 
Haematopinus suis [7]; (c) Culicoides melleus [14]; (d) Tenebrio molitor [11]; (e) Drosophila obscura [17]; (f) Tenebrio molitor [13]; (g) 
Aedes notoscriptus [15]; (h) Bacillus rossius [12]. All images reproduced with permission; (a) from S. K. Curtis and D. B. Benner, J. 
MorphoL, 210, 85-99, 1991, copyright 1991 John Wiley and Sons; (b) from B. Baccetti, Adv. Insect Physiol, 9, 315-397, 1972, copyright 
1972 Elsevier; (c) from J. R. Linley, Entomol. Exp. Appi, 26, 85-96, 1979, copyright 1979 John Wiley and Sons; (d) from B. Baccetti, 
A. G. Burrini, R. Dallai, F. Giusti, M. Mazzini, T. Renieri, F. Rosati, and G. Selmi, J. Mechanochem. Cell Motil, 2, 149-161, 1973, 
copyright 1973 Plenum Publishing Corporation, with kind permission from Springer Science+Business Media B.V; (e) from C. Bressac, 
D. Joly, J. Devaux, C. Serres, D. Feneus, and D. Lachaise, Cell Motil. CytoskeL, 19, 269-274, 1991, copyright 1991 John Wiley and Sons; 
(f) from D. M. Philips, in M. A. Sleigh, editor, Cilia and Fiagella, 379-402, 1974, copyright 1974 John Wiley and Sons; (g) from M. A. 
Swan, Gamete Res., 4, 241-250, 1981, copyright 1981 John Wiley and Sons; (h) from B. Baccetti, A. G. Burrini, R. Dallai, V. Pallini, P. 
Periti, F. Piantelli, F. Rosati, and G. Selmi, J. Ultrastruct. Res., 44, 1-73, 1973, copyright 1973 Elsevier. 



mistaken for an active, propagating wave under the mi- 
croscope. It has therefore been suggested that the sperm 
motility is caused solely by the minor wave. The relative 
extent of the contribution of the major and minor waves 
to propulsion is thus still an open question [14] . With the 
hydrodynamic study presented in this paper, we hope to 
provide physical insights on these unresolved problems. 

The structure of this paper is as follows. We idealize 
the double- wave structure as the propagation of superhe- 



lical waves and model the hydrodynamics using non-local 
slender body theory in §11. In §111, we present the com- 
puted hydrodynamic performance of spermatozoa of dif- 
ferent species and compare the predictions with available 
experimental data (§111 A). The features of super helical 
swimming arc illustrated by a specific model organism, 
namely the spermatozoa of Culicoides melleus (§111 B). 
We then investigate the effects of kinematic and geomet- 
ric parameters on the propulsion performance of a super- 
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helical swimmer (§111 C). Finally, the limitations of the 
present study and directions for future work arc discussed 
in §IV. 

II. MATERIALS AND METHODS 

A. Idealized double-wave structure: superhelical 
swimmers 

The experimentally observed double-wave structure 
of insect spermatozoa is mathematically idealized in 
this paper as a supcrhclix (a small helix itself coiled 
into a larger helix); we refer to the helical structure 
with the larger wavelength as the major helix and 
the other as the minor helix. To mathematically de- 
scribe a superhelix, we first construct the position vec- 
tor of a regular axial helix (the major helix) to be 
H(s') = [Am cos(kMCxs'), Am sui(kMtts'),as'], with a = 
1/ yJl + A 2 M k M . Here Um is the wave-number, Am is the 
amplitude, and s' <G [0,L'] and L' are the arc-length pa- 
rameter and the length of the major helix, respectively. 
From this basic helix the local unit tangent , unit nor- 
mal n^, and unit binormal = x vectors are 
determined, and we take them to define a local coordi- 
nate system upon which the minor helix is constructed 
[20]. The position vector of the combined superhelical 
shape is then given by 

X(s') = [x(s'),y(s'),z(s')] (1) 
= H(s') + A m cos(k m s')h A (s') 

±A m sm(k m s')b A (s'), (2) 

where A m and k m are the amplitude and wave-number 
of the minor helix, respectively. Two different configura- 
tions will be considered: the '+' sign leads to a superhe- 
lical structure where the major and minor helices both 
have the same chirality, whereas the ' — ' sign represents 
the case of opposite chirality. Note that s' is no longer 
the natural arc-length parameter, but merely a regular 
parameter for describing the swimmer's geometry. The 
arc-length of the complete superhelix, denoted by s, as 
a function of the parameter s' is determined by numeri- 
cal integration, and the total length of the superhelix is 
denoted by L. 

The major and minor helices are free to propagate 
waves at different wave speeds. Denoting cm and c m 



as the major and minor wave speeds respectively, the po- 
sition vector at time t, X(s',t), may be written in com- 
ponent form as 

x(s',t) = Am cos[fcA/(as' — CMt)} 

- A m cos[fc m (s' - c m t)] cos[k M (as' - c M t)] 
± aA m sin[fc m (s' - c m t)] sm[k M (as' - c M t)}, 

(3) 

y(s',t) = A M sm[k M (as' - c M t)\ 

- A m cos[fc m (s' - Cmt)] sin[k M (as' - c M t)] 

T ctAm sin[fc m (s' - c m t)] cos[k M (as - c M t)], 

(4) 

z(s', t) = ± aA m A M k M sin[fc m (s' - c m t)\ + as'. (5) 

In addition to the long flagcllum, insect spermatozoa 
have cell bodies which are very slender and short when 
compared to the flagellum size. We expect that the hy- 
drodynamic influence of the sperm head is negligible and 
do not include such a body in our consideration. A typ- 
ical superhelical swimmer is shown in Fig. 4, using the 
dimensionless parameters of Culicoides melleus sperma- 
tozoa. 

B. Hydrodynamical modeling 

Flagellar swimming is a result of the interaction be- 
tween the actuating body and the surrounding fluid. A 
tractable and accurate approach to studying such hy- 
drodynamic interactions exploits the slenderness of the 
flagellum, in which the velocity along the flagellar cen- 
terline is related to the fluid forces along the same curve. 
In previous studies, the fluid-body interaction has been 
modeled using a resistive force theory [21-23], in which 
local forces acting on the flagellum at any station along 
the filament are expressed in terms of the local velocity at 
the same location. Resistive force theory takes only local 
effects into account and neglects any hydrodynamic in- 
teractions between different parts of the deforming body. 
The local theory works well for simple geometries. How- 
ever, due to the complexity of the superhelical structure 
in the problem under consideration here, the local theory 
is inadequate (see Section HID for details), and we em- 
ploy instead the full non-local slender body theory [24] to 
study the hydrodynamics of the superhelical swimmers. 
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FIG. 4: Idealization of the double-wave structure as supcrhcliccs for the (a) same-chirality, and (b) oppositc-chirality configurations, using 
the dimensionless parameters of Culicoides melleus spermatozoa (Table II). 



The non-local theory captures the hydrodynamic inter- 
actions between distant parts of a curved filament, while 
still taking advantage of the slenderness of flagellum to 
simplify the analysis. 

The flagellum is modeled as a slender filament of length 
L and circular cross-section of radius e Lr(s), where e <C 
1 is the small aspect ratio of the flagellum (the maximum 
radius along the flagellum a r divided by its total length 
L), and r(s) is the dimensionless radius. The non-local 
slender body theory is algebraically accurate in ratio of 
the slenderness e; namely, by setting the radius profile 
to r(s) = y/4s(L — s)/L, the computed fluid velocity is 
accurate to 0(e 2 ) (see [24]). 

For a given velocity distribution v(s, t) along the fil- 
ament at time t, the corresponding fluid force per unit 
length f(s,i) is given implicitly by the non-local relation 



87r/!v(s,t) = -A[f(s,t)] -K[f(s,i)], 



(6) 



where 



A[f](s, t) = [c (I + tt) + 2(1 - tt)] f (s, t), (7) 
K[f](s,t) = (I + tt) t f ( g 'f)- f f 3 '*) d; 




o |s-s| 
I + RR I + ii 



R(s, s; t)\ \s — s\ 



f (s, t) ds, 



(8) 



are the local and non-local operators, respectively, fi is 
the shear viscosity of the fluid, cq = — ln(e 2 c) > 0, 
R(s,s;i) = X(s,i) - X(s,i), R = R/|R|, and i is the 
local unit tangent vector at the point s. 

The swimmers of interest deform their shapes in a 
prescribed, time-varying fashion (the superhelical wave 
pattern, Eqs. 3 to 5), and the velocity created on its 
surface by this deformation is given by Vd e form(s, t) = 
9X(s,t) / 'dt. At every time instant t, the swimmer can 
be seen as a solid body with unknown translational ve- 
locity U(t) and rotation rate Sl(t). The velocity cre- 
ated on the swimmer's surface due to swimming is then 
v sw im(s,£) = U + f2 x [X(s,t) — Xn], where X is an 
arbitrary reference point (taken as the origin here for 
simplicity). Therefore, the local velocity relative to the 
fluid v(s,i) is given by the sum of the deformation 
and swimming velocities: v(s, t) — Vdoform + v SW j m = 
U + n x [X(s,i) - X ] + dX(s,t)/dt. In this work, 
the wave propagation is towards the positive z direc- 
tion. Therefore, a negative swimming velocity (U) means 
the propulsion occurs in a direction opposite to the wave 
propagation, while a positive swimming velocity means 
both the propulsion and wave propagation occur in the 
same direction. 

Using a Galerkin method [25], we express the local 
force f (s, t) as a finite sum of Legcndre polynomials and 
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solve Eq. 6 for f(s) by requiring the equation to hold 
under inner products against the same basis functions. 
The first integral in the non-local operator K[f] is diag- 
onalized in this space [26, 27]. The system is closed by 
requiring the entire swimmer to be force free and torque 
free, 

/ f(s)ds = 0, (9) 
Jo 

I [X(s) - X ] x f (s)ds = 0, (10) 
Jo 

providing at each moment in time a system of 6 equations 
to solve for the 6 unknowns U(£) and H(t). 

The swimming velocities determined in the manner de- 
scribed above represent velocities in a reference frame 
fixed on the swimmer. In order to study the full 
three-dimensional swimming kinematics in the labora- 
tory frame (in which the body moves with velocity U(t) 
and rotation rate we must include a transformation 

between the two. We denote the Cartesian coordinate 
system moving with the swimmer, the body frame, as 
[e^e^ez] and the Cartesian coordinate system in the 
laboratory frame as [ei,e 2 ,e 3 ]. The evolution of the 
body frame, with respect to the laboratory frame, is then 
governed by 



dE 
~dt 



fl(t) x E, 



(11) 



where E = [e x ,e y ,e z ] T , along with the initial condition 
[e x ,e y ,e z ](t = 0) = [ei,e 2 ,e 3 ]. 



C. Non-dimensionalization 

The process of non-dimensionalization is very useful in 
science and engineering to identify the relevant dimen- 
sionlcss parameters governing the physics of the problem. 
In theoretical studies, it always allows a more concise de- 
scription of the system, while in experimental studies, it 
reduces the number of independent experiments required 
to fully explore the problem. The present system is made 
dimensionless by scaling lengths by l/k m , velocities by 
c m , and time by their ratio, l/c m k m . The dimensionless 
position vector describing the kinematics of the superhe- 



lix is thus given by 

x = R{cos[K(as' — ct)] — r cos(s' — t) cos[A'(as' — ct)] 
+ ar sin(s' - t) sm[K(as' - ct)]}, (12) 

y = R{sm{K(as' — ct)] — rcos(s' — t)sm[K(as' — ct)] 
— ar sin(s' — t) cos[K(as' — ct)]}, (13) 

z = as' + arR 2 K sin(s' - t), (14) 

where all variables are now understood to be dimension- 
less. Four dimensionless parameters characterizing the 
kinematics arc identified above: the dimensionless am- 
plitude of the major helix, R = AMk m , the ratio of the 
wave-numbers characterizing the major and minor he- 
lices, K = km /k m , the ratio of the minor helix ampli- 
tude to the major helix amplitude, r — A m j Am, and the 
ratio of the major wave speed to the minor wave speed 
c = Cm I Cm- In the double- wave pattern observed in in- 
sect spermatozoa, the major wave amplitude is always 
larger than the minor wave amplitude (r < 1; though 
this need not be true for a general superhelix). In ad- 
dition, it is also observed that the minor wave speed is 
always greater than the major wave speed in the double- 
wave structure of insect spermatozoa (c < 1). 

D. Kinematic and geometric data 

From the non-dimensionalization above, we have iden- 
tified four dimensionless parameters (r, R, K, c) required 
to fully characterize the centerline motion of a superhe- 
lical flagellum. Werner and Simmons [5] have compiled 
a very useful table containing kinematic and geometric 
data of the double- wave structure observed in insect sper- 
matozoa in previous studies. The table reveals that ex- 
perimental measurements of the necessary quantities for 
hydrodynamic modeling are very limited, due to the dif- 
ficulties involved in interpreting three-dimensional data 
from two-dimensional images [14]. 

Here we follow the table compiled by Werner and Sim- 
mons and estimate the missing information based on im- 
ages of insect spermatozoa reported in the literature. 
Table I contains reported and estimated data on the 
double- wave structure, including the wavelengths, ampli- 
tudes, frequencies, major and minor wave speeds, flagel- 
lum thickness, and the swimming speed of different in- 
sect spermatozoa. Our estimated quantities are marked 
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with stars to distinguish them from reported quantities 
by the original papers. In some studies (e.g. for Lygaeus 
[13]), only relative lengths can be given due to the lack of 
scale bars in the reported images; these relative quanti- 
ties are square-bracketed in Tables I & II. Geometric and 
kinematic data can display large variations even within 
a species {Tenebrio molitor [11]). When the distribution 
of the quantities are not given, arithmetic means of the 
available measurements are used whenever appropriate 
in the present study. 

The measurements are presented in the correspond- 
ing dimensionlcss quantities in Table II, using the scal- 
ings defined in §11 C. We have not found in the ref- 
erenced literature any information about the chirality 
of the major and minor helical structures. We there- 
fore present results below for both the same-chirality 
and opposite-chirality configurations. In our measure- 
ments (collected in Table I), the wavelength of the mi- 
nor wave is taken to be the two-dimensional distance be- 
tween adjacent minor wave peaks. A small correction 
factor is required to convert these two-dimensional quan- 
tities to the appropriate wavelengths in describing the 
three-dimensional supcrhclical structure. The correction 
factor depends on whether the supcrhclical structure is 
in the same-chirality (1/(1 — a\ m / Xm)) or the opposite- 
chirality (1/(1 + a\ m / \m)) configuration. 



III. RESULTS 



In this section, we first use the framework described 
above to predict the swimming performance of spermato- 
zoa of different species, and compare our hydrodynamic 
results with available experimental measurements. We 
then illustrate the basic features of superhelical swim- 
ming by focusing on a model organism, namely Culicoides 
melleus, and proceed to perform a parametric study in- 
vestigating the effects of certain kinematic parameters 
(minor and major wave speeds) and geometric parame- 
ters (minor and major wave amplitudes). 
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A. Hydrodynamic performance 
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Propulsion speed 

The propulsion speed is an important functional pa- 
rameter characterizing the motility of a sperm cell. There 
exist only very few measurements of the swimming speed 
of insect spermatozoa exhibiting the double-wave struc- 
ture: Ceratitis capitata was observed to swim with speed 
^sperm = 16 /im/s [16] and the minor wave speed is es- 
timated to be 150 /jm/s (hence, V — V^p C rm/cm=0.11); 
Culicoides melleus [14] was observed to swim with speed 
8.3 /jm/s and minor wave speed 80 /zm/s (V = 0.10). 
For the case of Tenebrio molitor, a wide range of sperm 
speeds were reported (from 16 /im/s to 100 /im/s) as 
well as minor wave speeds (40 /im/s to 300 /im/s) [11]. 
The distributions were not reported however, so using the 
arithmetic means we obtain a swimming speed to minor 
wave speed ratio of V = 0.34. 

For a superhelical swimmer in our model, the swim- 
ming kinematics are three-dimensional and unsteady in 
time [44], and the most relevant quantification of the 
propulsion speed is an average swimming velocity in the 
laboratory frame, denoted by U±, where the + and — 
signs represent the same- and opposite-chirality configu- 
rations respectively. The average swimming speed is de- 
fined as U± = \[{U X ), (U y ), (U z )]\, where ... denotes the 
magnitude of a vector, and (...) denotes a time average 
[45] . Since the chirality configuration remains unknown, 
we present predictions for both cases in Fig. 5a. Only the 
three sets of experimental measurements of the swimming 
performance (V = V spcrm /c rn ) mentioned above are avail- 
able for comparison and are super-imposed on the same 
figure (see Fig. 5 caption). Predictions as ratios of the 
swimming speed relative to the minor wave speed (di- 
mensionless speed) for species with no measurement of 
V are also provided in Fig. 5a. 

In most of the cases considered, the propulsion in 
the longitudinal direction (z-direction) of these super- 
helical swimmers is opposite to the direction of wave 
propagation (the major wave propagates in the positive 
z-dircction, and the minor wave propagates along the 
curved major helix and distally towards the positive z- 
direction), regardless of the chirality configuration. This 
is not unlike the behavior of swimmers propagating a 
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(a) U ± (b) V ± [%] 




FIG. 5: Predicted swimming performance: (a) average swimming speed, U±, of different species; (b) hydrodynamic efficiency, r)±, of 
different species. Open symbols (blue circles) represent the same-chirality configuration (+); filled symbols (red squares and diamonds) 
represent the opposite-chirality configuration (— ); for blue open circles and red squares, swimming occurs in the opposite direction as 
the wave propagation; on the contrary, red diamonds represent the cases of opposite-chirality configuration where the velocity in the 
2-direction occurring in the same direction as the wave propagation; black crosses represent experimental measurements of V = Vsperm/cm 
(See Tables I and fl). 



planar sinusoidal or a regular helical wave, for which the 
swimming direction is also opposite the direction of the 
wave propagation. In superhelical swimming, however, 
we also find cases of the opposite-chirality configuration 
where the body swims in the same direction as the wave 
propagation. Specifically, we note a qualitative difference 
between the study of the same-chirality and opposite- 
chirality configurations for Ceratitis capitata and Tene- 
brio molitor. In these simulations, a superhelical wave 
is set to propagate in the positive ^-direction, and both 
the opposite-chirality configurations of Ceratitis capitata 
and Tenebrio molitor generate a positive U z (swimming 
is in the same direction as the wave propagation), while 
their corresponding same-chirality configurations gener- 
ate a negative U z (swimming is in the opposite direc- 
tion as the wave propagation). The speed and efficiency 
of these peculiar swimmers are distinguished from other 
cases by red diamonds in Fig. 5. It remains a question 
whether this phenomenon may be observed in nature, 
since the chirality configuration and the swimming di- 



rection (relative to the wave propagation) of actual in- 
sect spermatozoa displaying the double-wave structure 
are still unclear. Nevertheless, this phenomenon by itself 
is intriguing and will be further explored in §111 C 1. 

The swimming speed predictions lie at least within 
the same order of magnitude of the experimental mea- 
surements for both chirality configurations. The same- 
chiralty results provide slightly better agreement than 
the opposite-chiralty results (without taking the swim- 
ming direction into account). For the same-chirality 
configuration, the discrepancies, \U+ — V\/V, between 
the predictions and the experimental measurements read 
36% for Ceratitis capitata, 49% for Culicoides melleus, 
and 18% for Tenebrio molitor. For the opposite-chirality 
configuration, the discrepancy between our predictions 
and the experimental measurements, |t/_ — V\/V, are 
65% for Ceratitis capitata, 50% for Culicoides melleus, 
and 21% for Tenebrio molitor. Given the primitive na- 
ture of the data employed (see §11 D), we consider the 
agreements here to be reasonable. However, we cannot 
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draw definite conclusions on the issue of chirality config- 
uration; further experimental observations are necessary. 

It shall be remarked that the present study is largely 
constrained by the unavailability of experimental mea- 
surement data. Critical kinematic information, such as 
the major and minor wave speeds, are often not reported 
in the literature and are impossible to estimate from the 
images available. The speed ratio c = cm / c m is not avail- 
able for simulations for most species. The computations 
here are still possible because of the independence of the 
swimming kinematics on the parameter c (verified nu- 
merically and will be explained in §111 CI). Therefore, 
the specific value of c is unimportant for all cases con- 
sidered here; we adopt c = in all simulations hereafter 
unless otherwise stated. 



2. Hydrodynamic efficiency 

Another important functional parameter is the hydro- 
dynamic efficiency of the swimmer. In the microscopic 
world, the hydrodynamic efficiency is typically very low. 
For a rigid helix, Lighthill [28] calculated theoretically the 
maximum efficiency attainable to be about 8.5%, while 
the typical efficiency of biological cells is around 1-2% 
[3, 28-32]. For low Reynolds number swimming, a com- 
mon measure of the hydrodynamic efficiency, ry, is the 
ratio of the rate of work required to drag the straight- 
ened flagellum through the fluid to the rate of work done 
on the fluid by the flagellum during swimming [28], 




where £|| = Att/cq is the drag coefficient for a straight, 
slender rod, and the brackets indicate a time average. 
We calculate the hydrodynamic efficiencies for both the 
same- and opposite-chirality configurations (see Fig. 5b) 
for spermatozoa of different species. For the same- 
chirality configuration, the efficiency ranges from 0.16% 
(Megaselia scalaris) to 1.2% (Aedes notoscriptus and 
Tenebrio molitor); for the opposite-chirality configura- 
tion, the efficiency ranges from 0.036% (Drosophila ob- 
scura) to 1.5% [Culicoides melleus). The efficiencies of 
these swimmers are comparable to typical biological cells. 



B. Model organism: Culicoides melleus 

In this section, we illustrate the features of supcrhelical 
swimming by singling out a supcrhelical swimmer defined 
using the geometric data of the sperm cell of Culicoides 
melleus (Table II). See Fig. 4 for the swimmer geometry. 

1. Swimming kinematics 

The swimming velocities computed for motion in the 
body frame [U x , U v , U z ] and in the laboratory frame 
[U x ,U y ,U z ] of the superhelical swimmer with data from 
the sperm cell of Culicoides melleus are plotted in Fig. 6 
for the same-chirality (Fig. 6a) and opposite-chirality 
(Fig. 6b) configurations. We include in these figures the 
results for ratio of the major and minor wave speeds c = 
(solid lines) and c = 0.4 (dotted lines). When c = 0, 
the deformation (Eqs. 12 to 14) is 27r-periodic, there- 
fore the swimming kinematics in the body frame are also 
27r-periodic. However, when observed in the laboratory 
frame, the coupling between the translational and ro- 
tational kinematics renders the swimming velocities no 
longer 27r-pcriodic, and the motion is unsteady in time 
(Fig. 6c & d). We observe a pattern consisting of higher 
frequency oscillations modulated by a lower frequency 
envelope. For the case of c = 0.4, when observed in 
the body frame, we see modulated waveforms. However, 
when transformed to the laboratory frame, the cases of 
c = and c = 0.4 have identical swimming kinematics 
(see the overlapped solid and dotted lines in Fig. 6c & 
d), implying that the swimming kinematics in the labo- 
ratory frame are independent of the major propagating 
wave speed. The value of c affects only the kinematics in 
the body frame; why the relative wavespeed is unimpor- 
tant is described in greater detail in §111 C 1. 

The three-dimensional swimming velocities give rise to 
a doubly-helicated trajectory (the presence of a minor 
structure on top of a major helical structure). How- 
ever, this could be difficult to observe experimentally, 
since the major amplitude of the doubly-helicated tra- 
jectory is usually much smaller than that of the supcr- 
helical swimmer; the swimmer would apparently move 
with a straight trajectory (with very small oscillations in 
the transverse direction). Recall that for regular helical 
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FIG. 6: Three-dimensional swimming velocities in the body frame U = (U x , U y , U z ) for (a) the same-chirality and (b) the opposite-chirality 
configurations, and in the laboratory frame U = (U x ,Uy,Uz) for (c) the same-chirality and (d) the opposite-chirality configurations, for 
Culicoides melleus spermatozoa. 



swimming (r — 0), the trajectory of the helical swimmer 
reduces to a regular helix. 

2. Head-less swimming 

The shape and size of the sperm head vary among sper- 
matozoa of different species: human and bull spermato- 
zoa have relatively large, paddle-shaped heads, whereas 
insect spermatozoa have elongated heads which are al- 
most indistinguishable from the mid-piece. The ad- 
ditional hydrodynamic resistance from the presence of 
a head would seemingly degrade the swimming perfor- 
mance of the sperm cell. However, Chwang and Wu [33] 
showed that a sperm head is actually necessary for helical 
swimming; without one, the motion is that of a rotating 
rigid body, which cannot be realized absent an external 



force or torque [34]. To satisfy the zero net-torque con- 
dition, a sperm head is required to balance the reaction 
torque acting on the fiagellum. This constraint does not 
apply for planar, sinusoidal wave motion, which can swim 
without an anchor or load. 

We pause to point out a subtle but important differ- 
ence between a rotating prokaryotic helical tail and a 
cukaryotic tail propagating a bending helical wave. For 
a cukaryotic tail propagating a bending wave, the fluid 
forces act to rotate the fiagellum opposite the direction 
of the apparent helical rotation. The rotation due to the 
fluid reaction creates torques due to local spinning (rota- 
tion of the fiagellum about its center line), which balance 
the opposing torque generated by the helical wave prop- 
agation. Therefore, a cukaryotic cell could theoretically 
swim without a sperm head, albeit very slowly because 
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the flagellum is very slender and the local torques are 
correspondingly small. Quantitatively, using a local drag 
model ([33]) it can be shown that the head- less swim- 
ming speed scales as U/c ~ 2fik 2 + k 2 A 2 )b 2 + 0(b 4 ), 
where b is the radius of a cross-section of the flagellum, 
and k and A are the wavenumber and helical radius re- 
spectively. Using the geometrical data of flagella of Eu- 
glena viridis summarized by Brennen and Winet [4] , and 
assuming a flagellar diameter of 2b « 0.25/im [35, 36], 
we find U/c « 1CP 3 . For a prokaryotic tail, however, 
the helix rotates as a rigid body. In this case, the local 
spinning torques generated by the active helical rotation 
and the passive rotation due to the fluid reaction are 
identical in magnitude but opposite in sign. The torque- 
free condition therefore requires that the fluid reaction 
counter-rotates the helix at precisely the rotation rate 
of the helical wave propagation. Hence there can be no 
effective helical rotation, and the body cannot swim. Al- 
though head-less swimming is theoretically possible for 
eukaryotic tails, the swimming speed would be excep- 
tionally small as shown by the estimation above. The 
subtle difference between the two types of helical waves 
just described is often therefore neglected, and it is gen- 
erally reasonable to state that head-less swimming is not 
possible using helical wave propagation. 

In superhelical swimming a sperm head is not required 
for self-propulsion (indeed, all cases reported in this pa- 
per were studied with the absence of a sperm head) since 
a superhelical wave motion is in general not a rigid body 
motion. Furthermore, since actual insect spermatozoa 
heads are very slender and short compared with the entire 
length of the flagellum, the contribution of its hydrody- 
namic resistance and the hydrodynamic interactions with 
the flagellum should be negligible. Therefore, we do not 
expect the presence of a slender and short sperm head to 
introduce qualitative differences in the results. 



C. Parametric study 

Next, we explore the effects of certain kinematic (wave 
speeds) and geometric parameters (wave amplitudes) on 
the swimming velocities of a superhelical swimmer. 



1. The effect of c 

In this problem, we have scaled the velocities upon 
the minor wave speed c m , which implies that the dimen- 
sionlcss swimming velocities scale linearly with the mi- 
nor wave speed. Here we examine the effect of the major 
wave speed on the swimming performance of a superhe- 
lical swimmer and answer the question: to what relative 
extent do the two waves contribute to the propulsion of 
the superhelical swimmer [14]? Specifically, we study the 
effect of the parameter c = cm/ c mi which is the ratio of 
the major wave speed to the minor wave speed. For every 
species shown in Table II, we fix all parameters but vary 
the value of c from zero to unity (c < 1, as cm < c m ). It 
is found that the swimming velocities in the laboratory 
frame is independent of the value of c. We have already 
shown in Fig. 6 the results for two values of c for illus- 
tration. While the resulting swimming motions in these 
two cases differ significantly in the body frame, they are 
identical when observed in a laboratory frame of refer- 
ence. The numerical results imply that the major wave 
speed does not contribution to propulsion. 

These findings may be understood by noting an ambi- 
guity in definition. First consider a single helical filament 
(without a sperm head) placed in a fluid. Such a body 
cannot swim on its own, for there is no load or cell to 
counter balance the torque it would exert on the sur- 
rounding fluid during rotation, so it must be motionless 
as seen in a fixed, laboratory frame. However, this body 
may be represented as an active helical body propagat- 
ing a wave with velocity c plus a rigid body rotation that 
contributes a wave with speed — c. That the body frame 
may be chosen arbitrarily allows for such an ambiguity 
in the definition of the swimming speed, but in the lab- 
oratory frame this ambiguity disappears. 

A similar argument can be used to show that c has 
no bearing on the swimming speed of a superhelical flag- 
ellum in the laboratory frame. For a superhelical flag- 
ellum, the propagation of the major wave can be de- 
fined as a rigid body rotation of the entire superhelical 
structure about the longitudinal axis (z-axis) in the body 
frame. The apparent rotation rate is the sum of the rota- 
tion caused by the active propagation of the major wave 
and the rotation caused by fluid reaction to maintain 
the force-free and torque-free conditions. For swimming 
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without a head, only the apparent rotation rate is impor- 
tant because the fluid forces and torques are functions of 
the apparent rotation rate alone (not the absolute value 
of the rotation rate due to the active major wave prop- 
agation). Any alternation of the major wave speed is 
accompanied by a corresponding change in the rotation 
rate caused by the fluid reaction, resulting in the same 
apparent rotation rate and the same force and torque bal- 
ances. Therefore, the kinematics and dynamics do not 
depend on the absolute value of the major wave speed c 
(or c m ). In other words, the major wave speed does not 
contribute to swimming. 

For the case of swimming with a head, the fluid reac- 
tion will not only rotate the flagcllum but also the sperm 
head, which creates extra torques and perturbs the orig- 
inal force and torque balances. Altering the absolute 
value of the major wave speed in this case theoretically 
will affect the swimming velocities even in the laboratory 
frame, because it changes the relative portion of the ro- 
tation rate caused by the fluid reaction and hence the 
value of the extra torque from the sperm head. However, 
since the sperm head is slender and short compared with 
the overall length of the flagellum of actual insect sper- 
matozoa, the extra resistant forces and torques created 
should be insignificant. Therefore, we suggest that the 
propagation of the major wave contributes very little to 
the propulsion (none in the case of head- less swimming) , 
and that it is the minor wavespeed which is primarily re- 
sponsible for the propulsion. Note also that the hydrody- 
namic efficiencies in §111 A 2 are independent of the value 
of c. Therefore, from a hydrodynamic efficiency point of 
view, there is no advantage or disadvantage to actively 
propagate a major wave. However, there are other energy 
costs not taken into account here (for example, work done 
to produce the sliding of microtubules within the flagel- 
lum) in actively propagating a major wave. Therefore, 
we speculate that it might be energetically more favor- 
able for a doubly-hclicated organism to propagate only a 
minor wave. 

2. The effect of R and r 

We now examine the effect of geometrical dimension- 
less parameters R and r, keeping all other parameters 



of Culicoides melleus spermatozoa fixed. R is the di- 
mensionless major wave amplitude and r is the ratio of 
minor to major wave amplitude, hence the minor wave 
amplitude is given by rR (see the schematic in Fig. 7). 
In general, one expects the propulsion speed to increase 
with the wave amplitude for simple geometries. However, 
for a supcrhclical structure, the dependence of the aver- 
age swimming speed on the minor wave amplitude dis- 
plays interesting behavior. Physically, both the propul- 
sive force and the bulkincss of the structure are varied 
upon changing the major or minor wave amplitudes. The 
competition between these factors and the coupling be- 
tween kinematics in different directions create interesting 
geometric dependencies of the average swimming speed. 
We illustrate this by observing the average swimming 
speed under different frames of reference. 

First, we look at the time-averaged swimming velocity 
in the body frame U h = \[{U x (t)), (U y (t)), (U z (t))]\. As 
shown in Fig. 7a(i), for the same-chirality configuration 
the average swimming velocity grows monotonically with 
R and r. Note that, keeping other parameters fixed, in- 
creasing the value of R = AMk m for a fixed r = A m j Am 
geometrically means that both the major Am and mi- 
nor A m amplitudes are increased simultaneously by the 
same proportion. There is a competition between the in- 
crease in the overall hydrodynamic resistance due to the 
increased bulkincss (correlated with increases in Am and 
A m ) and an enhanced propulsive force (correlated with 
an increase in A m ). In the case of the same-chirality con- 
figuration, the latter effect dominates. However, for the 
opposite-chirality configuration, as shown in Fig. 7a(ii), 
non-monotonic variations in the swimming speed with 
the major amplitude R are observed for certain values of 
r, the ratio of the minor to major wave amplitudes. 

Since the swimming kinematics are three-dimensional, 
variations of the geometric parameters affect swimming 
velocities and rotational rates in all directions. In partic- 
ular, a large propulsion speed in the body frame does not 
necessarily imply a large net propulsion speed in the lab- 
oratory frame. The mean swimming speeds in the labora- 
tory frame for both chirality configurations are shown in 
Fig. 7b(i-ii). The coupling between the swimming kine- 
matics in all directions produces more complicated vari- 
ations in the propulsion speed as a function of R and r. 
Non-monotonic behaviors are observed in both cases. 
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FIG. 7: Parametric study of the dependence of the average swimming speed in the body fixed frame U\, (a(i)&(ii)) and in the laboratory 
frame U± (b(i)&(ii)) as a function of the dimcnsionless parameters R and r (see the schematic for geometrical illustration). Panels c(i) 
& c(ii) show the average swimming velocity in the ^-direction in the body frame, (U z ), as a function of R and r. The panels on the left 
(right) refer to the same- (opposite-) chirality configuration. The dotted line in panel c(ii) represent the contour of (U z ) = 0. Geometric 
data of Culicoides melleus spermatozoa are used for other fixed parameters. 
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(a) (b) 




(i) (ii) (i) (ii) 

FIG. 8: Three-dimensional swimming velocities in the laboratory frame ((a)i & (b)i) of two opposite-chirality superhelical swimmers, and 
their corresponding geometries ((a)ii & (b)ii) . (a) R = 7, r = 0.4, and U z is negative, while the wave propagation is towards the positive 
^-direction; (b) R = 9, r = 0.4, and U z is positive, while the wave propagation is towards the positive 2-direction. Geometric data of 
Culicoides melleus spermatozoa arc used for other parameters. 



We have already noted that some superhelical swim- 
mers propel themselves surprisingly in the same direc- 
tion as the wave propagation, unlike for planar or sin- 
gle helical wave propulsion. For the range of param- 
eters explored in this paper, this direction reversal is 
found to occur only in the opposite-chirality configura- 
tion for sufficiently large R and r (Fig. 7c(ii)), while the 
average propagation velocity (U z ) is always negative for 
the same-chirality configuration (i.e. the swimming direc- 
tion is always opposite to the direction of the propagat- 
ing wave) (Fig. 7c(i)). In Fig. 8, we show the detailed 
swimming velocities of two opposite-chirality superheli- 
cal swimmers, where one of them has its longitudinal 
propulsion in the opposite direction relative to the wave 
propagation (R = 7, r = 0.4, Fig. 8a), and the other 
in the same direction as the wave propagation (i? = 9, 
r = 0.4, Fig. 8b). It is intriguing that a small change 
in the geometry (see the corresponding superhelices in 
Fig. 8a(ii) & b(ii)) can lead to a swimming direction re- 
versal. We argue in the next section (§111 D) that this 
transition is related to the hydrodynamic interaction be- 
tween distinct parts of the superhelical flagellum. 



D. Comparison between slender body theory and 
resistive force theory 

In order to consider the relative importance of non- 
local hydrodynamic interactions in the swimming of su- 
perhelices, we now compare our results to those obtained 
using the more commonly used local drag model (eq. (6), 
but neglecting the non-local term K). The local drag 
model (so-called resistive force theory) ignores hydrody- 
namic interactions between distinct parts of the curved 
flagellum, and is expected to work well [21-23] for simple 
geometries where different parts of the body are suffi- 
ciently well separated. However, for more complicated 
geometries, the local drag model may not capture even 
the correct qualitative features. Recently, Jung et al. 
[20] studied the rotational dynamics of opposite-chirality 
superhelices and found a qualitative discrepancy (the ro- 
tational direction of the superhelix when being towed in 
a viscous flow) between the experimental results and the 
predictions from the resistive force theory. Other recent 
works have also shown the inadequacy of the local drag 
model [32, 37]. 

Figure 9 shows the swimming speeds computed with 
both non-local slender body theory and resistive force 
theory, using two sets of geometrical data: Culi- 
coides melleus (Fig. 9a (same-chirality) & b (opposite- 
chirality)), and Tenebrio molitor (Fig. 9c (same- 
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FIG. 9: Comparison between the resistive force theory (RFT) (dashed lines) and the non-local slender body theory (SBT) (solid lines), 
using geometric data of Culicoides melleus spermatozoa (a) & (b), and Tenebrio molitor spermatozoa (c) &i (d). The panels on the left 
(right) refer to the same- (opposite-) chirality configuration. 



chirality) & d (opposite-chirality)). In these two cases, 
all other parameters are fixed but the amplitude ratio 
r is varied from to 0.5. Increasing r complicates the 
swimmer geometry, and the hydrodynamic interactions 
are expected to be more significant for large r. For the 
case of Culicoides melleus (Fig. 9 a & b), we see good 
(even quantitative) agreement between the results form 
the slender body theory (solid lines) and the resistive 
force theory (dashed lines). However, for the case of 
Tenebrio molitor (Fig. 9 c & d), which has larger values of 
R and K, the deviation between the two models becomes 
significant for the samc-chirality configuration (Fig. 9c) 
as r increases. There are even qualitative discrepancies: 



in the opposite-chirality configuration (Fig. 9c) for large 
values of r, the resistive force theory fails to capture the 
transition in the sign of (U z ) predicted by the slender 
body theory. Since no such transition is found when 
the hydrodynamic interactions are ignored, this transi- 
tion may be attributed to nonlocal hydrodynamic inter- 
actions of the body with itself. In general, and perhaps 
unsurprisingly, we have shown that the local drag model 
breaks down when the geometry of the structure is suffi- 
ciently intricate. 
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E. Asymptotic analysis 

In the locomotion of some species of spermatozoa, the 
minor wave amplitude is much smaller than the major 
wave amplitude (see Tables I & II), which motivates us 
to perform an asymptotic analysis for r <C 1. Such an 
asymptotic analysis linearizes the problem geometrically 
and allows the nonlinear effects to be taken into account 
order by order, making the problem more amenable to 
mathematical analysis. However, even in the asymp- 
totic consideration, three-dimensional force and torque 
balances do not yield tractable analytical results. Hence, 
in the spirit of Chwang and Wu [33] , we perform the force 
and torque balances only in the longitudinal (z-) direc- 
tion using the resistive force theory, which is expected to 
be at least qualitatively correct in the asymptotic limit 
r <C 1. The local force f z and torque m z may be ex- 
pressed as 

fz = c UzFz (t, s)U z + cn z F z (t, s)Cl z + c Fz (t, s), (16) 
m z = c UzMz (t,s)U z + cn zMz (t,s)Q z +c Mz (t,s), (17) 

where the coefficients arc determined analytically from 
Eq. 6 (without the non-local operator) and the geome- 
try of the swimmer (Eqs. 12 to 14). We consider regular 
perturbation expansions in r for every term in Eqs. 16 
and 17, and enforce the force-free and torque-free con- 
ditions order by order. A non-zero time-averaged swim- 
ming velocity enters at 0(r 2 ). The leading order mean 
swimming velocities for the cases of same- and oppositc- 
chirality configurations read 



U± 



R r 



1J2. 



{(2-2L 2 )(£-l) 



2L 2 (l + R 2 K 2 ) 3/2 
+ R 2 K 2 [2(£ - 1 ± K) - L 2 (2£ - 2 ± K)} 

-2 [£- l + R 2 K 2 {i- l±K)] cosLX, (18) 

where £ = £_i_/£|| is the ratio of the drag coefficients 
in the normal direction to the longitudinal direction 
(£|| = 47t/c ,£_l = 87r/(2 + c )). For very slender fila- 
ments, £ « 2. We note that the substitution K — > —K 
in Eq. 18 converts the same-chirality speed U+ to the 
opposite-chirality speed [/_ . The propulsion speed has a 
quadratic dependence on the minor wave amplitude for 
small r, which is also found in planar [38] and helical [39] 
geometries. 



These asymptotic results (Eq. 18) are compared to the 
finite amplitude simulations (both the resistive force the- 
ory and the slender body theory predictions) in Fig. 10, 
for a very small amplitude ratio of r = 0.01 and a slen- 
derness ratio of e = 1/1000. There is excellent agreement 
between the asymptotic results and the finite amplitude 
resistive force theory simulations. The discrepancy be- 
tween the non-local slender body theory and the local 
drag model highlights the importance of non-local hy- 
drodynamic interactions for such organisms. 

We conclude by pointing out an intriguing theoretical 
curiosity, that drag anisotropy (£ ^ 1) is not required 
for superhelical swimming: setting £ = 1 in Eq. 18, non- 
zero mean propulsion velocities (of equal magnitude but 
opposite signs) are obtained for both chirality configura- 
tions, 



- R*K 3 {L 2 -2 + 2cosL) 2 

±U ' T 2L(1 + R 2 K 2 fl 2 T ' 



(19) 



Note that Becker et aUs [40] argument of the require- 
ment of drag anisotropy for locomotion is only true for 
inextensible swimmers and does not apply here. The su- 
perhelical kinematics described (Eqs. 12 to 14) are only 
possible when extensibility is allowed; the minor helix 
is built upon another curved structure (the major helix) 
and local extension and contraction is implied in the wave 
kinematics. When extensibility is permitted, the relax- 
ation of the drag anisotropy requirement has been re- 
cently shown [41]. That the swimming speed is non-zero 
is due to intrinsic variations in length (and hence drag) 
embedded in the curved geometry of the superhelices. 
A minor helix built upon a major helix has relatively 
shorter lengths in the regions closer to the longitudinal 
axis, creating an overall imbalance of hydrodynamic drag 
even in the isotropic drag case (£ = 1). A similar exam- 
ple is a toroidal helix (a helix built upon a circle) , which 
is an idealized model studied recently for dinoflagcllates 
[42] . We expect that the propagation of a wave along a 
toroidal helix should also require extensibility, and that 
propulsion is still possible even without drag anisotropy 
[411. 
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FIG. 10: Comparison of the asymptotic results (dotted lines) with the finite-amplitude simulations using the resistive force theory (dashed 
lines) and the non-local slender body theory (solid lines). Shown are the contour lines of the average swimming velocity U± for the (a) 
same-chirality, and (b) opposite-chirality configuration for a small minor to major wave amplitudes ratio of r = 0.01. 



IV. DISCUSSION 

In this paper we have studied a morphologically inter- 
esting double-wave structure exhibited by various insect 
spermatozoa. The construction of such spermatozoa is 
considerably more complex than those for flagella which 
exhibit simpler planar or helical waves: the flagellum 
does not only have a more complicated 9+9+2 arrange- 
ment of microtubules but also mitochondrial derivatives 
and accessory bodies running along the axoneme [5] . We 
have mathematically idealized the double- wave structure 
as a superhelical structure and presented a hydrodynamic 
study on superhelical swimming. The available data is 
primitive and sparse; nevertheless, we consider the agree- 
ment between experimental measurements and the the- 
ory explored herein to be quite reasonable. Through 
numerical experiments, we have found that the major 
wave speed has little contribution to propulsion when the 
sperm head is small, as is the case for insect spermatozoa. 
When there is no sperm head, the propulsion speed is in- 
dependent of the major wave speed and depends entirely 
upon the minor wave speed. We have also explored the 
dependence of the propulsion speed on the dimensionlcss 
major wave amplitude R and the ratio of the minor to 
major wave amplitudes r (Fig. 7), and counter-intuitive 



behaviors have been found for the opposite-chirality con- 
figuration. In particular, we have found that propulsion 
and wave propagation can occur in the same direction for 
supcrhcliccs in the opposite-chirality configuration. 

The present study suggests that the major wave has 
negligible influence on the motility of a superhelical 
swimmer. This finding favors the recent hypothesis by 
Werner et al. [10] that the major helical wave is a static 
(non-propagating) structure; the minor wave structure is 
solely responsible for the motility, and the apparent ma- 
jor wave propagation is simply due to the passive rotation 
of the entire geometry (see §1). However, in the study by 
Baccetti et al. [16] on the motility of Ceratitis capitata, 
they have adopted the same experimental techniques as 
in Gibbons et al. [43] , which distinguished the rolling fre- 
quency from the apparent beat frequency. In their work, 
the rolling frequency was measured by stroboscopic ob- 
servation of the eccentrically attached sperm head, and 
the flagellar beat frequency was measured by the same 
means with the sperm head adhered to the bottom of 
the observation dish. Using this method, the major wave 
speed measured should be taken as an active propaga- 
tion speed. We are not in a position to provide a definite 
answer on whether or not the major wave propagates 
actively in actual insect spermatozoa. However, accord- 
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ing to the present study, we suspect that there might be 
some biological reasons, other than motility, for the ma- 
jor wave to propagate actively. We also do not know if 
the propagation of the major wave is a biological pre- 
requisite for the propagation of the minor wave. Further 
biological studies are required to answer these questions. 

It is illuminating to compare the superhelical swim- 
ming studied here to regular helical swimming. The 
swimming trajectories arc qualitatively different in the 
two cases: in regular helical swimming, the trajectory 
is a regular helix, whereas in superhelical swimming the 
trajectory is doubly-helicated. In addition, a regular he- 
lical flagcllum cannot swim on its own; a sperm head is 
required to swim (since the absence of a sperm head ren- 
ders the deformation of the regular helix a rigid body 
motion). In contrast, the propagation of superhelical 
waves along a flagellum is in general not equivalent to 
rigid body motion, and hence "head-less" swimming is 
possible. This might help explain the presence of only 
a small slender sperm head in insect spermatozoa. In 
other words, superhelical swimming can be viewed as an 
alternative mechanism to regular helical swimming when 
only a small sperm head is available. 



Finally, a non-local slender body theory was used in 
this work and compared with a simpler and widely used 
local drag model. We showed that the resistive force 
theory failed to capture even the qualitative features of 
the swimmer when the geometry becomes complicated. 
The results suggest that further hydrodynamic studies 
on superhelical structures require more advanced models 
than the local resistive force theory, such as the slender 
body theory employed here. Optimization with respect 
to the efficiency of the swimmer taking into account of 
the viscous dissipation, and other energy costs due to 
bending and internal sliding of filaments [31], which have 
been neglected in this work, will be interesting for future 
work. 
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